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Abstract 

Vigna reflexo-pilosa, which includes a neglected crop, is the only one tetraploid species in genus Vigna. The ancestral species 
that make up this allotetraploid species have not conclusively been identified, although previous studies suggested that a 
donor genome of V. reflexo-pilosa is V. trinervia. In this study, 1,429 azul<i bean EST-SSR markers were developed of which 38 
EST-SSR primer pairs that amplified one product in diploid species and two discrete products in tetraploid species were 
selected to analyze 268 accessions from eight taxa of seven Asian Vigna species including V. reflexo-pilosa var. glabra, V. 
reflexo-pilosa var. reflexo-pilosa, V. exilis, V. hirtella, V. minima, V. radiata var. sublobata, V. tenuicaulis and V. trinervia to 
identify genome donor of V. reflexo-pilosa. Since both diploid and tetraploid species were analyzed and each SSR primer pair 
detected two loci in the tetraploid species, we separated genomes of the tetraploid species into two different diploid types, 
viz. A and B. In total, 445 alleles were detected by 38 EST-SSR markers. The highest gene diversity was observed in V. hirtella. 
By assigning the discrete PCR products of V. reflexo-pilosa into two distinguished genomes, we were able to identify the two 
genome donor parents of Creole bean. Phylogenetic and principal coordinate analyses suggested that V. hirtella is a species 
complex and may be composed of at least three distinct taxa. Both analyses also clearly demonstrated that V. trinervia and 
one taxon of V. hirtella are the genome donors of V. reflexo-pilosa. Gene diversity indicates that the evolution rate of EST- 
SSRs on genome B of Creole bean might be faster than that on genome A. Species relationship among the Vigna species in 
relation to genetic data, morphology and geographical distribution are presented. 



Citation: Chankaew S, Isemura T, Isobe S, Kaga A, Tomooka N, et al. (2014} Detection of Genome Donor Species of Neglected Tetraploid Crop Vigna reflexo-pilosa 
(Creole Bean), and Genetic Structure of Diploid Species Based on Newly Developed EST-SSR Markers from Azuki Bean {Vigna angularis). PLoS ONE 9(8}: el 04990. 
doi:l 0.1 371/journal.pone.Ol 04990 

Editor: Manoj Prasad, National Institute of Plant Genome Research, India 
Received April 29, 2014; Accepted July 16, 2014; Publislied August 25, 2014 

Copyriglit: © 2014 Chankaew et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All azuki bean EST-SSR primers, along with the 
corresponding SSR motif, product size, and primer sequence, are available at http://marker.kazusa.or.jp/Azuki. All other relevant data are within the paper and its 
Supporting Information files. 

Funding: The work was supported by the Thailand Research Fund (TRF). The funders had no role in study design, data collection and analysis, decision to 
publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* Email: tomooka@affrc.go.jp (NT}; agrpss@yahoo.com (PS) 
9 These authors contributed equally to this work. 



Introduction 

The Leguminosae genus Vigna comprises about 100 species. 
These species are morphologically diverse and geographically 
widespread and mainly found in Africa (African Vigna; subgenus 
Vigna) and Asia (Asian Vigna; subgenus Ceratotropis) [1]. The 
Asian Vigna comprises 2 1 species in which seven are domesticated 
and/or cultivated in various geographical and climatic regions, 
and cropping systems in Asia [2]. The seven domesticated/ 
cultivated species include V. aconitifolia (moth bean), V. angularis 
(azuki bean), V. mungo (black gram), V. radiata (mungbean), V. 
reflexo-pilosa (creole bean), V. stipulacea (jungli bean) and V. 



umhellata (rice bean). All Asian Vigna species are diploid having 
1 1 haploid cliromosomes (2n = 2x = 22) with the exception for 
Creole bean which is a tetraploid species with number of haploid 
chromosome of 22 (2n = 4x = 44). In fact, Creole bean is the only 
natural amphidiploid in the subtribe Phaseolinae [3]. Cytogenetic 
analyses of V. reflexo-pilosa showed that the species formed 22 
bivalents without multivalent at meiosis and was, therefor 
considered to be an amphidiploid [4] [5]. 

Cultivated and wild forms of Creole bean are classified as V. 
reflexo-pilosa var. glabra and V. reflexo-pilosa var. reflexo-pilosa, 
respectively [6]. Wild creole bean is widely distributed in East, 
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Southeast and South Asia, and across the islands from the west to 
the north Pacific islands. It is also found in Papua New Guinea and 
northern Austraha [7] [8] [9]. The cultivated Creole bean was 
formerly recognized as a glabrous variety of mungbean, V. radiala 
var. glabra [1]. Then, it was treated as a distinct species, V. 
glabrescens [3]. It is differentiated from its wild progenitor 
principally by thick glabrous stem and erect growth habit and 
reported to be cultivated as pulse in Vietnam and the Philippines 
or as forage in India, Mauritius and Tanzania [2] . This crop shows 
resistance to several insect pests and diseases such as bruchids, 
bean fly, powdery mildew, and cucumber mosaic virus [10], and is 
partially cross-compatible with mungbean [1 1] . Thus, Creole bean 
has potential to be a gene source for breeding other Vigna crops. 
In addition, cultivated Creole bean could be considered as a novel 
crop for the future and wild Creole bean as a wider genepool to 
improve cultivated Creole bean. 

The origin and genome donors of Creole bean have been the 
subject of debate. Based on the study on isozyme banding patterns 
and interspecific hybridization among the Ceratotropis species, 
Egawa et al [10] [12] proposed that V. reflexo-pilosa var. reflexo- 
pilosa is the descendant of natural interspecific hybridization 
between V. trinervia and V. miriima followed by spontaneous 
chromosome doubling. An acc(\ssion used as V. minima in Egawa 
et al. [10] [12] was correctly identified as V. hirtella by Konarev et 
al. [13], and was used in the present study (No. 21, JP108851, 
from Malaysia). Proteinase inhibitors polymorphism study in Asian 
Vigna hy Konarev et al. [13] supported that one genome donor of 
Creole bean is V. trinervia and suggested that the other genome 
donor is most likely V. hirtella or its closely related species. The 
phylogenetic studies based on plastid sequence data also supported 
that V. trinervia is a genome donor of Creole bean [14] [15]. In 
contrast, rDNA-ITS sequence variation showed high similarity 
between Creole bean and V. exilis, V. hirtella and V. umhellata 
[16], suggesting that one of these species is the genome donor of K. 
reflexo-pilosa. 

The objective of this study was to determine the putative 
genome donor species of tetraploid Creole bean using EST-SSR 
markers. To do so, we developed EST-SSR from azuki bean and 
used them to analyze accessions of Creole bean and other Vigna 
species that are candidate genome donor species. 

Materials and Methods 

Plant Materials 

A Japanese azuki bean cultivar 'Erimo-shouzu' iV. angularis 
var. angularis, accession no. JP377.^2) obtained from the Gene- 
bank, National Institute of Agrobiological Sciences (NIAS), 
Tsukuba, Japan was used for development of EST-SSR markers. 
Eight accessions consisting of four major Vigna crop species; azuki 
bean, rice bean, black gram and mungbean were used for assessing 
transferability and polymorphisms of the EST-SSR markers 
(Table 1). Two hundred and sixty-eight accessions from eight 
taxa of seven Asian Vigna species (Figure 1, Table SI) including 7 
of V. reflexo-pilosa var. glabra, 51 of F. reflexo-pilosa var. reflexo- 
pilosa, 13 of F. exilis, 47 of V. hirtella, 49 of V. minima, 13 of F. 
radiata var. sublohata, 42 of V. tenuicaulis, and 46 of V. trinervia 
were used to analyze the genome origin of V. reflexo-pilosa. 

Development of Azuki Bean EST-SSR Markers 

Total RNA was extracted from seedlings and young pods of 
Erimo-shouzu using Plant RNA Purification Reagent (Invitrogen, 

CA, USA). Purification of polyadenylated RNA and conversion to 
cDNA were performed as described by [1 7] . Synthesized cDNA 
was resolved by 1% agarose gel electrophoresis, and fragments 



ranging from 1 to 3 kb were recovered. The recovered fragments 
were cloned into the Eco 'Sl-Xho I site of the pBluescript II SK- 
plasmid vector (Stratagene, CA, USA) and introduced into the E. 
coli ElectroTen-Blue strain (Stratagene, CA, USA) by electropo- 
ration. To generate ESTs, plasmid DNAs were amplified from the 
colonies using TempliPhi (GE Healthcare UK Ltd, Buckingham- 
shire, England) and subjected to sequencing using a BigDye 
Terminator Cycle Sequencing Ready Reaction Kit (Applied 
Biosystems, CA, USA). The reaction mixtures were run on an 
automated ABI PRISM 3730 DNA Analyzer (Applied Biosys- 
tems). 

Sequencing chromatograms were converted into nucleotide 
bases with Phred [18] [19], and the sequences derived from the 
vector and linkers were removed with CROSSMATCH [19]. The 
EST reads were quality-trimmed with TRIM2 [20] using the 
Phred quality score S20, and ambiguous regions including more 
than ten X or N bases were trimmed. Contiguous, high-quality 
reads > 100 bp were submitted to the DDBJ/EMBL/GenBank 
databases under the accession numbers HX939204 to HX950377 
(11,174 entries). The PHRAP program with default parameters 
was used to cluster and identify non-redundant azuki bean ESTs 
[19]. 

Simple sequence repeats (SSRs) Sl5 nucleotides in length, 
which contained all possible combinations of di-nucleotide (NN), 
tri-nucleotide (NNN), and tetra-nucleotide (NNNN) repeats, were 
identified from the non-redundant azuki bean ESTs using fuzznuc 
program in EMBOSS [21] for SSRs within two mismatches. 
Primer pairs for amplification of SSR-containing regions were 
designed based on the flanking sequences of each SSR with the aid 
of the PrimerS [22] so that the amplified fragment sizes were 
between 90 bp and 300 bp in length. The newly developed 
markers were designated as VES (Vigna EST-derived SSR) 
markers. 

EST-SSR Marker Analysis 

Total genomic DNA was extracted from young leaves of each 
Vigna accession using the method described by [23] with a slight 
modification. The DNA was quantified against a lambda DNA on 
1.0% agarose gel stained with ethidium bromide and diluted to 
5 ng/|j.l for PGR amplification. 

Amplification, transferability and polymorphisms of all the 
EST-SSR markers were confirmed using eight accessions, two 
each from azuki bean complex, blackgram, rice bean and 
mungbean (Table 1). Five microliters of PGR reaction mixture 
including 0.5 ng of total genomic DNA, 0.02 U BIOTAQDNA 
polymerase (BIOLINE, UK), Ix PGR buffer (BIOLINE, UK), 
3 mM MgCla, 0.2 mM dNTPs and 2 pmol of the forward and 
reverse EST-SSR primers. The PGR thermal cycling was 
performed as follows: 94°G for 1 min; 3 cycles of 94°C for 30 s 
and 70°G for 30 s, followed by 3 rounds of the same program in 
which the annealing and extension temperatures were decreased 
by 2°G every 3 cycles; 3 cycles of 94°G for 30 s, 62°G for 30 s and 
72°G for 30 s, followed by 2 rounds of the same program in which 
the annealing temperatures were decreased by 2°G every 2 cycles; 
30 cycles of 94°G for 30 s, 55°G for 30 s and 72°C for 30 s, and a 
final cycle at 72°C for 10 min. The PGR products were separated 
by 10% polyacrylamide gel electrophoresis in tris-borate-ethylene 
diamine tetraacetic acid (TBE) buffer according to the standard 
protocol, and banding patterns in the gel were stored as pictures. 
The characteristics of PGR products such as intensity, banding 
pattern and size range were recorded for each accession and EST- 
SSR markers. Polymorphism between wild and cultivated 
accessions in each species was recorded. 
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^ Vigna exilis 

^ V. reflexo-pilosa var. glabra 

^ V.tiirtella 

I V. minima 

■ V. reflexo-pilosa var. reflexo-pilosa 



I V. radiate var. sublobata 
^ V. tenuicaulis 
^ V. trinervia 



Figure 1. Distribution of 286 accessions of 8 taxa of tKie genus Vigna subgenus Ceratotropis. Nunnbers of accessions analyzed are shown 
within circles, triangles and squares. 
doi:1 0.1 371 /journal.pone.01 04990.g001 



Out of 1,429 azuki bean EST-SSR markers, 175 markers with 
good transferability in all the four major Vigna crop species were 
labeled with fluorescent dyes and further examined for polymor- 
phism in a panel of eight taxa of Asian Vigna consisting of 1 6 
accessions, two accessions for each species randomly chosen from 
the 268 accessions (Table SI). The 5 '-end of the reverse primers 
were labeled with one of the following fluorescent dyes; 6-FAM 
(blue), HEX (green), and NED (yellow) (Applied Biosystems). Five 
microliters of PGR reaction mixture contained 5 ng of genomic 
DNA, 1 X QIAGEN Multiplex PGR Master Mix, and 5 pmol of 
the forward and reverse primers. PGR reactions were performed 
in a GeneAmp PGR System 9700 (Apphed Biosystems). The PGR 
thermal cycling was programmed as foUows: 95°G for 15 min 
followed by 40 cycles of 94°G for 30 s, 55°G for 90 s, 72°G for 
60 s, and a final cycle at 72°C for 10 min. One microliter of 10 
times dilution PGR product was mixed with 8.5 |XL of Hi-Di 
formamide and 0.125 |J,L of ROX size standard (Applied 
Biosystems). The mixture was denatured at 95°G for 5 min and 
then run on an ABI PRISM 3130x1 DNA Analyzer (Applied 
Biosystems). Allele size for the highest stutter peak with the height 
ranging between 500 and 10,000 RFU were recorded. The 
genotyping was performed using GeneMapper 4.0 (Applied 
Biosystems). 

Thirty-eight EST-SSR primer pairs amplified one product in 
diploid species and two discrete products in tetraploid species with 
height ranging RFU (Table S2) were selected for further analysis 



in the 268 Vigna accessions. Three primers with different labels 
and size products were mixed into a single PGR reaction mixture 
and amplified as multiplex PGR. After genotyping, a single allele 
size was scored for each marker in each accession corresponding to 
the strongest peak. 

Data Analysis 

Since both diploid (1 genome set) and tetraploid (2 genome set) 
species were analyzed together in this study and each SSR primer 
pair detected 2 loci in the tetraploid species, we separated the 
genome of V. reflexo-pilosa, into 2 different diploid types, viz. A 
and B. Therefore 10 genomes were recognized from the 8 taxa of 
7 Vigna species. 

Genetic distance (Z)a) [24] for all possible pairs of accessions was 
calculated using software POPULATIONS 1.2.28 (available at 
www.cnrs-giffr/pge/bioinfo/populations). Da among the 10 
genomes was also calculated using the same software. A 
phylogenic tree was constructed to reveal relationships among 
accessions based on Da using neighbor-joining clustering method 
by software MEGA 5.05 [25] with bootstrap support (1,000 
rephcates) obtained by re-sampling the allelic frequency data. In 
order to confirm the results from neighbor-joining clustering, 
principal coordinate analysis (PGoA) was also performed using 
software PAST [26] to reveal the relationship among different 
accessions. In addition, based on sources of V. hirtella, Seehalak et 
al. (2006) divided V. hirtella in to 3 subgroups (al, a2 and b), thus 
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the relationship of individuals based on PCoA analysis enabled 12 
sub-genome groups to be determined. Using average Da genetic 
distance between 12 sub-genomes, an unrooted dendrogram 
showing relationships between these genomes was constructed by 
neighbor-joining method using software POPULATIONS 1.2.32. 
(available at www.cnrs-giffr/pge/bioinfo /populations) with boot- 
strap support (400 replicates) obtained by re-sampling the allelic 
frequency data. The trees were visualized with MEGA ver. 5.05 
[25]. 

Results 

Azuki Bean EST-SSR Markers 

A total of 11,167 cDNA clones were sequenced consisting of 
7,534 clones from a seedling library and 3,633 clones from a 
young pod library. After clustering, 4,896 potential non-redundant 
EST sequences, including 2,350 contigs and 2,546 singletons, were 
generated with a total of 4,284,693 qualified bases (Table S3). By 
using the fuzznuc program in EMBOSS, a total of 1,188 SSRs 
were identified in the 4,896 non-redundant EST sequences. Of the 
1,188 SSRs, di-, tri-, and tetra-nucleotide SSRs accounted for 
71.6%, 26.7%, and 1.7%, respectively (Table S4). Assuming that 
total length of the non-redundant azuki bean EST sequences is 
4.3 Mbp, the frequency of occurrence of the SSRs in transcribed 
regions was estimated to be one in every 3.6 kb. Altogetherl05 
primer pairs were initially designed on the flanking regions of the 
132 perfect SSR motifs (Table S4). To increase the number of 
candidate EST-SSR markers, additional 1,324 primer pairs were 
designed on the flanking regions of 1,545 imperfect SSR motifs 
allowing one or two base mismatching. As a result, a total of 1,429 
EST-SSR markers were designed (Table S4), of which 149, 28, 5 
and 5 markers identified 2, 3, 4 and £5 SSRs in the regions 
between the primer parrs, respectively. Thus the total number of 
identified SSRs by the 1,429 EST-SSR markers was 1,677 that 
consisted of 137 (8.2%) di-nucleotide repeats, 1,400 (83.5%) tri- 
nucleotide repeats and 140 (8.3%) tetra-nucleotide repeats (Table 
S4). Among the di-nucleotide repeats, poly(AG),-i (n = 88, 5.2% of 
total) were most frequently observed, followed by poly(AT),-, 
(n = 35, 2. 1 %) and poly (AC)„ (n = 14, 0.8%). Among the ten types 
of tri-nucleotide repeats observed, poly(AAG)„ (n = 322, 19.2%) 
were the most abundant, followed by poly(GGA)n (n = 195, 1 1.6%) 
and poly(ATC)„ (n=180, 10.7%). Among the thirteen tetra- 
nucleotide repeats, poly(AAAG),-, (ii = 45, 2.7%) were the most 



frequently observed, followed by poly(AAAT)„ (n = 30, 1.8%), and 
poly(AAAC)„ (n = 20, 1.2%). The details of the designed azuki 
bean EST-SSR primers, along with the corresponding SSR motif, 
product size, and primer sequence, are available at http://marker. 
kazusa.or.jp/Azuki and in Table S5. 

Amplification and Transferability of Azuki Bean EST-SSR 
Markers 

Transferability and polymorphism of the 1,429 azuki bean 
EST-SSR markers were initially examined by polyacrylamide gel 
electrophoresis using eight accessions consisting of four major 
Vigna crop species; azuki bean complex, black gram, rice bean 
and mungbean (Table 1, Table SI). These wild and cultivated 
parental accessions had been used for a linkage map construction 
in each species ([27] [28] [29] [30]. Amplification of azuki bean 
EST-SSR markers in the same species, V. angularis, and the 
closely related species, V. umhellata, were 91.3 and 92.9% which 
were higher than those in V. radiata and V. mungo (Table 1). 
However, the size and pattern of PGR products of EST-SSR 
markers did not always reveal simple banding pattern with 
expected size, even in the same species. Among the amplified 
markers, 81.8 to 83.6% of them revealed simple banding pattern 
suitable for further application. When polymorphism between 
cultivated and wild accessions was examined, V. radiata possessed 
the highest polymorphism in which 277 markers (23.1% of 
amplified markers) were polymorphic. Further, 1 75 markers with 
good transferability in the four Vigna species were tested for 
polymorphism in a panel of 1 6 accessions from 8 taxa of Asian 
Vigna (Table SI). The test revealed that 22 markers (12.6%) failed 
to amplify, 36 (20.6%i) amplified some accessions, whereas 117 
(66.8%) successfully amplified all the 16 accessions (Table 2). 
Among the 1 1 7 amplifiable markers, 2 (1.1%) gave product size 
larger than 500 bp, 4 (2.3%) amplified multiple products, 45 
(25.7%) amplified one product in tetraploid species, 10 (5.7%) 
were monomorphic, and 56 (32%) were polymorphic, amplifying 
one product in diploid species and two discrete products in 
tetraploid species and were considered to be suitable for analyzing 
origin of genome donor of the tetraploid species. 

EST-SSR Polymorphism and Genetic Diversity 

Thirty-eight primer pairs of azuki bean EST-SSR were 
analyzed in 268 accessions of 10 genomes from eight taxa of 
Asian Vigna. Based on discrimination of the diploid species we 



Table 2. Characteristics of 196 EST-SSR primers developed in this study. 




Type 


Description 


Number of primers (%) 


1 


Not amplified in all accessions 


22 (12.6) 




Amplified in some accessions 


36 (20.6) 


2-1) 


1) Amplified in 1 to 8 accession(s) 


5 (2.9) 


2-2) 


2) Amplified in 9 to 15 accessions 


31 (17.7) 




Amplified in all 16 accessions 


117 (66.8) 


3-1) 


1) PCR product of more than 500 bp 


2 (1.1) 


3-2) 


2) IVlultiple PCR product 


4 (2.3) 


3-3) 


3) Monomorphic PCR product 


10 (5.7) 


3-4) 


4) Single PCR product in tetraploid 


45 (25.7) 


3-5) 


5) Two PCR products in tetraploid 


56 (32) 


Total 




175 (100.0) 


doi:1 0.1 371 /journal.pone.Ol 04990.t002 
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distinguished two discrete products detected by eacli marker of V. 
reflexo-pilosa into two types, viz, A and B. The product size that is 
the same or very similar to V. trinervia was scored as type A, while 
the other product size was scored as type B. In total, 445 alleles 
were detected in 10 genomes by the 38 EST-SSR loci (Tables 3, 
4). The number of alleles detected per locus ranged between 3 
(VES0777 and VES1271) and 34 (VESl 172) with a mean of 1 1.7 
alleles per locus (Table 4). The PIC values ranged from 0.36 
(VES0116) to 0.94 (VESl 172) with a mean of 0.67. None of the 
EST-SSR markers had PIC values lower than 0.3, while 20 
markers had PIC values higher than 0.7. The markers showed 
high PIC value in V. hirtella (0.44) and V. minima (0.42), but low 
value in V. reflexo-pilosa (0.04-0.09). In general, each marker 
showed a PIC value of 0 in both populations of V. reflexo-pilosa. 
The allelic richness was between 2.0 for marker VES0777 and 9.9 
for marker VESl 172 with a mean of 4.7. Fifteen primers had 
allelic richness higher than 5.0 (Table 4). 

V. hirtella possessed the highest gene diversity (0.478), while V. 
reflexo-pilosa showed very low gene diversity (0.053 to 0.094) 
(Table 3). All genomes of V. reflexo-pilosa showed no observed 
heterozygosity, while the other genomes showed very low observed 
heterozygosity from 0.008 to 0.028 (Table 3). 

Genetic Relationship among Genomes 

Genetic distance (-Da) among the ten genomes is shown in 
Table 5a. In most cases Dj^ among genomes was high (>0.6). 
Nonetheless, Da between V. trinervia and V. reflexo-pilosa var. 
glabra (A) or V. reflexo-pilosa var. reflexo-pilosa (A) was low being 
0.333 and 0.313, respectively. While, both V. reflexo-pilosa var. 
glabra (B) and V. reflexo-pilosa var. reflexo-pilosa (B) showed 
lowest Da with V. hirtella being 0.585 and 0.589 (Table 5a) or 
with V. hirtella (al), being 0.429 and 0.422 (Table 6), respectively. 
Interestingly, V. trinervia, which belongs to section Angulares 
showed lower Da with V. radiata var. sublobata (section 
Ceratotropis) than with those species in the section Angulares. 

A neighbor-joining tree was constructed based on the genetic 
distances for all possible pairs of 268 accessions. The tree showed 
two major groups with 96.9% bootstrap value, namely azuki bean 
group and mungbean group (Figure 2). V. exilis, V. minima, V. 
tenuicaulis, V. hirtella, V. reflexo-pilosa var. glabra (B), and V. 
reflexo-pilosa var. reflexo-pilosa (B) were clustered in the azuki 
bean group. While V. radiata var. sublobata, V. trinervia, V. 



reflexo-pilosa var. glabra (A) and V. reflexo-pilosa var. reflexo- 
pilosa (A) were clustered in the mungbean group. In the azuki 
bean group, 13 accessions of V. exilis and 49 accessions of V. 
minima were independent clusters with 95.3% and 83.5% 
bootstrap values, respectively. V. hirtella showed high divergence 
in which 9 accessions in the population of V. hirtella (b) clustered 
strongly with V. tenuicaulis with 91.5% bootstrap value, while the 
other 38 accessions in V. hirtella (a) clustered with V. reflexo-pilosa 
var. glabra (B) and V. reflexo-pilosa var. reflexo-pilosa (B) with 
83.8% bootstrap value. Those 38 accessions of V. hirtella also 
showed two subclusters al and a2 (Figure 2 and 3). In the 
mungbean group, V. trinervia clustered with V. reflexo-pilosa var. 
glabra (A) and V. reflexo-pilosa var. reflexo-pilosa (A) with a 100% 
bootstrap support. Five accessions from Myanmar (40-44) formed 
a distinct branch with a bootstrap value of 87.8% and were 
separated from the other accessions (Figure 2). Three accessions 
from Laos (31-33) were differentiated from the other accessions 
and were clustered together with V. reflexo-pilosa (A) accessions. 

The principle coordinate analysis (PCoA) was also conducted to 
confirm species relationship based on the neighbor-joining tree 
using the same genetic distance estimates. The first three PCs 
together accounted up to 56.6% of the total variation. The first, 
second and third PCs accounted for 38.1%, 10.6% and 7.9%, 
respectively (Figure 4). A three-dimensional PC plot (PCI, PC2 
and PC 3) of the twelve genome groups is shown in Figure 4, which 
clearly separated the genomes into 5 distinct groups. Group I 
comprised V. reflexo-pilosa var. reflexo-pilosa (A), V. reflexo-pilosa 
var. glabra (A), and V. trinervia. Group II comprised V. minima 
only. Group III comprised V. reflexo-pilosa var. glabra (B), V. 
reflexo-pilosa var. reflexo-pilosa (B), V. exilis and V. hirtella (al, 
a2). Group FV comprised V. hirtella (b), and V. tenuicaulis. Group 
V comprised solely V. radiata var. .sublobata. Distribution range of 
V. reflexo-pilosa var. glabra (A) and V. reflexo-pilosa var. reflexo- 
pilosa (A) in Group I was narrower than V. reflexo-pilosa var. 
glabra (B) and V. reflexo-pilosa var. reflexo-pilosa (B) in Group III. 
Of all species used in this study, only V. hirtella was unambig- 
uously distinguished into two sub-groups, a and b. Subgroup a 
showed close relationship with populations of V. reflexo-pilosa var. 
reflexo-pilosa (B); V. reflexo-pilosa var. glabra (B) and V. exilis, 
while subgroup b showed close relationship with V. tenuicaulis 
(Figure 3). 



Table 3. Genome number assigned In this study with number of alleles, gene diversity and observed heterozygosity analyzed by 
38 EST-SSR markers. 



Genome 


No. of accessions 


No. of loci typed 


No. of alleles 


Gene diversity 


Observed heterozygosity 


1 [Vigna exilis) 


13 


38 


100 


0.305 


0.012 


2 [V. reflexo-pilosa var. glabra (a)] 


7 


38 


44 


0.053 


0.000 


3 [V. reflexo-pilosa var. glabra (b)] 


7 


38 


44 


0.052 


0.000 


4 ty. hirtella) 


47 


38 


171 


0.478 


0.019 


5 {V. minima) 


49 


38 


203 


0.458 


0.017 


6 [V. reflexo-pilosa var. reflexo-pilosa (a)] 


51 


38 


56 


0.058 


0.000 


7 [V. reflexo-pilosa var. reflexo-pilosa (b)] 


51 


38 


76 


0.094 


0.000 


8 (V. radiata var. sublobata) 


13 


38 


101 


0.357 


0.008 


9 {V. tenuicaulis) 


42 


38 


122 


0.309 


0.022 


10 (I/, trinervia) 


46 


38 


98 


0.306 


0.028 


Total 


326 


38 


445 


0.711 


0.013 


doi:1 0.1 371 /journal.pone.Ol 04990.t003 
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Discussion 

Transferability of Azuki Bean EST-SSRs 

One of the advantages of SSRs derived from EST (EST-SSR 
markers) to other DNA marker systems is the liigh transferability 
across species and genera because of primer pairs are designed 
from conserved transcribed regions, although they are generally 
less polymorphic than genomic SSR markers (gSSR). In this study, 
more than 83% of the azuki bean EST-SSR markers were able to 
amplify DNA from 8 Asian Vigna species including four major 
crop species, viz. V. angularis, V. umhellata, V. radiata and V. 
mungo of the subgenus Ceratotropis . This result is in agreement 
with that of [31] who reported that more than 80% of mungbean 
EST-SSR markers were transferable to other species in the 
subgenus Ceratotropis. Transferability of the EST-SSR markers 
from azuki bean (in this study) and mungbean [31] is greater than 
that of gSSR markers from azuki bean (>67% transferable; [28]) 
and mungbean (>59% transferable; [32]). The high transferability 
of azuki bean EST-SSRs and of mungbean EST-SSRs to the 
related Vigna species indicates high genome conservation among 
the species in the subgenus Ceratotropis . Therefore, the azuki bean 
EST-SSR markers developed in this study will be useful for 
comparative genomic and genetic diversity studies in Vigna 
species. Nonetheless, as expected, polymorphism rate of the azuki 
bean EST-SSR markers detected between cultivated and wild 
accessions of azuki bean, rice bean, mungbean and black gram is 
low (14% in black gram and 23% in mungbean). In the same plant 
materials, Chaitieng et al. [28] reported that polymorphism rate of 
the genomic azuki bean SSRs was as high as 50%) in black gram to 
63% in mungbean. Thus low polymorphism rate of the azuki bean 
EST-SSRs can be an undesirable character for the use of these 
markers in comparative genome mapping in the genus Vigna. 
However, EST-SSRs represent true genetic diversity which may 
be more directly associated with traits of interest in breeding as 
compared to gSSRs [33], and they may reflect better relationships 
among the related species in genetic diversity study. 

Evolution and Domestication of Tetraploid V. reflexo- 
pilosa 

V. reflexo-pilosa is the only tetraploid species of the genus 
Vigna. Although previous studies on Vigna species using DNA 
markers and plastid sequences clearly demonstrated that V. 
trinervia is a genome donor of V. reflexo-pilosa, those studies 
were unable to identify the other genome donor. In the present 
study, five putative genome donor species of V. reflexo-pilosa 
belong to the section Angulares and one species representing the 
section Ceratotropis were analyzed with azuki bean EST-SSRs. By 
assigning the discrete PGR products of V. reflexo-pilosa into two 
distinguished genomes (designated as genome A and B), we were 
able to identify the two genome donor parents of Creole bean. 
Phylogenetic tree (Figure 2) and PCoA plot (Figure 4) clearly 
demonstrated that genome A of V. reflexo-pilosa is received from 
V. trinervia. This result confirms the previous results reported by 
[7] [10] [13] [14] [15] that V. trinervia is a diploid genome donor. 
V. trinervia and V. reflexo-pilosa shares several similar morpho- 
logical characters such as seed shape (rectangular), non-protruded 
hUum and large golden flower. The phylogenetic tree and the 
PCoA plot also unambiguously demonstrated that the other 
genome (genome B) of V. reflexo-pilosa is originated from V. 
hirtella. This result is in agreement with the result reported by 
Tateishi [7] [10] [13] [16] that V. hirtella is one of candidate 
genome donor of Creole bean. Results from phylogenetic analysis 
of plastid DNA sequences of [14] [15] suggested that V. trinervia 
is the maternal genome donor of V. reflexo-pilosa. Since the 
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Figure 2. A phylogenetic tree showing relationship among 286 accessions in 12 sub-genome groups of the genus Vigna subgenus 
Ceratotropis based on variation at 38 EST-SSR loci. 

doi:10.1371/journal.pone.0104990.g002 



previous and our results demonstrated close genetic relationship 
between tetraploid V. reflexo-pilosa and diploid V. trinervia/V. 
hirtella, we propose that Creole bean evolved from interspecific 
hybridization between V. trinervia as female parent and V. hij'tella 
as male parent, followed by genome duphcation. 

Phylogenetic tree (Figure 2) and PCoA plot (Figure 4) showed 
lower divergence of genome A of V. reflexo-pilosa compared to 
genome B. Very low gene diversity within both genomes A and B 
suggested that this tetraploid species has a monophyletic (single) 
origin and only evolved recently. Gene diversity within the 
genome A in both V. reflexo-pilosa var. reflexo-pilosa (wild form) 
and V. reflexo-pilosa var. glabra (cultivated form) is very similar 
(0.058 vs. 0.053), while gene diversity of the genome B in wild 
Creole bean is about twice of that in the cultivated one (0.94 vs. 
0.52) (Table 3). This indicates that the evolution rate of EST-SSRs 
on genome B of Creole bean might be faster than that on genome 
A. The marked morphological differences between the cultivated 
and wild Creole beans are thicker and erect glabrous stem in the 
former. Moreover, cultivated V. reflexo-pilosa still possesses a 
relatively high degree of pod shattering (Somta and Chankaew, 
personal observation), an important domestication trait for legume 
crops [34]. Therefore, the cultivated V. reflexo-pilosa can be 
treated as a semi-domesticated form. Based on cultivation, 



utilization and distribution of diploid genome donor species, 
Creole bean appears to have been domesticated in Southeast Asia. 
The crop is now very rarely cultivated, although it is produced in 
northern mountainous villages in Vietnam under the same name 
and consumed in the same way as mungbean [2]. 

Genetic Structure of Diploid Asian Vigna Species 
Detected by EST-SSR 

We found that V. hirtella possessed greater genetic variation 
than V. exilis, V. minima, and V. tenuicaulis. A previous study 
based on AFLP suggested that V. hirtella germplasm collected 
from Thailand, Malaysia and Myammar is a species complex 
consisting of two taxa (types), called V. hirtella (a) and V. hirtella 
(b) [35]. Although the study used only 3 accessions for V. hirtella 

(a) and 5 accessions for V. hirtella (b), aU the 3 V. hirtella (a) 
accessions were included in V. hirtella (al) group, and 5 V. hirtella 

(b) accessions were included in V. hirtella (a2) group in the present 
study. In addition, our EST-SSR results detected the existence of 
another genetically distinct type of materials designated as V. 
hirtella (b) which is closely related to V. tenuicaulis (Figure 2, 3, 4). 
These accessions should be re-examined morphologically whether 
they can be included within a morphological variation of V. 
tenuicaulis. 
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Figure 3. A phylogenetic tree showing relationship among 12 sub-genome groups of the genus Vigna subgenus Ceratotropis 
constructed based on the genetic distance shown in Table 6. 

doi:1 0.1 371 /journal.pone.01 04990.g003 



V. trinervia is widely distributed across Asia and also found in 
Papua New Guinea, Madagascar and East Africa. It is the second 
most widely distributed species in the subgenus Ceratotropis after 
V. radiata var. sublobata [2]. Nevertheless, gene diversity of V. 
trinervia was not high as compared to the other related species 
(Table 3). In Southeast Asia, especially in Thailand and the 
peninsular Malaysia, V. trinervia is often found on rural roadside 
habitats and in or near to rubber or oil palm plantations. It seems 
that some of the populations in those areas were introduced 



recently as a cover plant in the plantations [2] . This may account 
for the low diversity of V. trinervia. 

V. tenuicaulis distributes in open wet habitats in northern 
Southeast Asia, and appeared to be closely related to V. hirtella (b) 
and V. angularis [13] [35] [6] [36]. Tomooka et al. suggested that 
there are no major barriers to hybridization between V. 
tenuicaulis and V. hirtella (b) [2] . Our results confirmed the close 
relationship between these two species (Figure 2, 4). Previous 
studies based on molecular, biochemical and morphological 
variations all showed high level of distinctness and variation 
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Figure 4. PCA scattered plot depicting relationship among 286 accessions in 12 sub-genome groups of the genus Vigna subgenus 
Ceratotropis based on variation at 38 EST-SSR loci. 

doi:1 0.1 371 /journal.pone.01 04990.g004 



within V. tenuicaulis. EST-SSR marker variation in our study 
revealed contrasting results. Although we used as many as 30 
accessions of V. tenuicaulis from three countries, viz. Laos, 
Myanmar and Thailand. (Figure 1, Table SI), the species showed 
relatively low gene diversity (0.309). 

V. exilis has been reported only in Thailand and Myanmar [2] 
[14]. In Thailand, this species is restricted to rocky limestone 
mountains. Its habitats suggest that it may be useful as gene 
sources for resistant to alkaline soil and drought conditions. SSR 
analysis revealed high level of intra-specilic diversity in V. exilis 
[37]. Our results also supported this despite only 13 accessions 
from narrow geographical origin (west of Thailand) of this species 
were used (Figure 1, Table SI), the species showed similar level of 
gene diversity to V. tenuicaulis and V. trinervia (Table 3). 

V. minima has broad environments adaptation and grows well 
in shaded deciduous forest floors and open-wet habitats in East 
and Southeast Asia [2] . It is the only species in section Angulares 
that is found on the forest floor. Among the Asian Vigna species 
analyzed in this study, V. minim.a is the second most diverse species 
after V. hirtella (Table 3). The high difiFerentiation of V. minima is 
due to wide geographical distribution of the analyzed accessions. 
In Southeast Asia, isolation of V. minima in forests in different 
mountainous regions across Thailand and Myanmar and its 
sporadic occurrence in patches of forests in those regions may 
account for high level of population divergence [35]. Accessions of 



V. minima in the present study were most closely related with V. 
hirtella, especially V. hirtella (b), followed by V. exilis (Figure 2, 4). 
Similar finding was reported by [35]. This agrees with their 
morphological appearance [2] that V. minima is sometimes 
confused with V. hirtella. V. minima can be distinguished from 
V. hirtella by smaller bracteole and more protruding hUum with 
weU-developed rim-aril [2]. 

In summary, one of the advantages of EST-SSR markers to 
other DNA marker systems is their co-dominant nature and high 
transferability to closely related species. In this study, we used 
EST-SSR primer pairs, which amplified one PGR product in 
diploid species and two discrete PGR products in tetraploid species 
to separate genomes of V. 1'eflexo-pilosa into 2 different diploid 
types, viz. A and B, that could successfully determine dilferent sets 
of the donor genomes. Phylogenetic analyses revealed that V. 
trinervia and V. hirtella are the genome donors of V. reflexo- 
pilosa. Both genomes of cultivated and wild Creole bean accessions 
showed low divergence suggesting that domestication of V. reflexo- 
pilosa is a relatively recent event. 

Supporting Information 

Table SI Origin of 268 accessions used to analyze origin 
of genome of V. reflexo-pilosa. 

(XLSX) 
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Table S2 Transferability and polymorphisms of azuki 
bean EST-SSR markers in the other Vigna species. 

(XLSX) 

Table S3 EST composition in non-redundant azuki 
bean EST sequences. 

(XLSX) 

Table S4 Summary of SSRs in the non-redundant azuki 
bean ESTs and designed EST-SSR primers. 

(XLSX) 



Table S5 Primer sequences, SSR motifs, expected sizes 
and EST sequences of the designed EST-SSR markers. 

(XLSX) 
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